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Abstract. In the present paper, an integrable semi-discrete analogue of the one-dimensional 
coupled Yajima-Oikawa system, which is comprised of multicomponent short-wave and one 
component long-wave, is proposed by using Hirota’s bilinear method. Based on the reductions 
of the Backlund transformations of the semi-discrete BKP hierarchy, both the bright and dark 
soliton (for the short-wave components) solutions in terms of pfaffians are constructed. 
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1. Introduction 

Both the nonlinear Sehrodinger (NLS) equation and the Yajima-Oikawa (YO) equation are 
import models arising in a variety of physieal eontexts lUl . It is known that they ean be derived 
from the so-ealled k-eonstrained KP hierarehy [|2l. Speeifieally, the NLS equation eorresponds 
to the k-eonstrained KP hierarehy with k = I, while the YO equation is assoeiated with the 
ease with k = l [|3l. 

The integrable diseretization of nonlinear Sehrodinger equation 


^^n,l — (1 “h 0|(5'n I ) (^n+1 T 1) 


( 1 ) 


was originally derived by Ablowitz and Ladik and was named Ablowatiz-Ladik (AL) 

lattiee. Similar to the eontinuous ease, the AL lattiee equation admits bright soliton solution 
for the foeusing ease (o = 1) [IHISl and dark soliton solution for the defoeusing ease (o = — 1) 
Il6l . The geometrie eonstruetion of the AL lattiee was given by Doliwa and Santini |I71. 

The semi-diserete eoupled nonlinear Sehrodinger equation 



( 2 ) 


(3) 


where O; = ±1(/ = 1,2), is of importanee both mathematieally and physieally. It was solve 
by the inverse seattering method in dHl HI. The general multi-soliton solution in terms of 
pfaffians was found reeently in ifTOlfTTI . whieh is of bright type for the foeusing-foeusing ease 
(Oi = 02 = 1), is of dark type for the defoeusing-defoeusing ease (Oi =02 — —1), and eould 
be of mixed type for the foeusing-defoeusing ease (Oi = 1, 02 = — 1). 

In eompared with the sueeessful eonstruetion of the integrable diserete analogue of the 
NLS equation and the eoupled NLS equations, the integrable diserete analogues of the YO 
equation and its multi-eomponent generalizations are missing. It should be pointed out 
that fully diserete NLS and YO equations were eonstrueted most reeently ifT^ . however, 
their semi-diserete limits eannot eonverge to the model sueh as the Ablowatiz-Ladik lattiee. 
Therefore, it is of important to eonstruet an integrable semi-diserete analogue of the one¬ 
dimensional (ID) eoupled Yajima-Oikawa (YO) system: 


=0, A' = 1,2, • • • ,M, 

M 


(4) 


(5) 


where are arbitrary real eonstants, and L indieate the /Jth short-wave (SW) and long 
wave (LW) eomponents, respeetively. 
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When M = 1, the original YO system deseribed sueh a physical process, in which a 
resonant interaction takes place between a weakly dispersive LW and a SW packet when 
the phase velocity of the former exactly or almost matches the group velocity of the latter 
[fT4]| . Yajima and Oikawa ifTSl proposed the model equation for the interaction of a Langmuir 
wave with an ion-acoustic wave in a plasma and showed that the system is integrable by 
the inverse scattering transform method. This model equation was also derived in diverse 
physical backgrounds such as hydrodynamics [IT^ . nonlinear optics [[171 UHl and biophysics 
[[T9]l . It admits both bright and dark soliton solutions [12011211 l22l . The rogue wave solutions 
to the ID YO system have recently been derived by using Hirota’s bilinear method ll23l 
and Darboux transformation [l2^ . For the multicomponent ID YO system ©-([Sl), it can 
be obtained from the KP-hierarchy by reduction as shown in Ref. [I251l2^l271l . Most recently, 
Kanna et al. [|28l considered the multicomponent coupled ID YO system , and pointed 
out that the coupled YO system with two component LWs (M = 2) can be deduced from 
a set of three coupled NLS equations governing the propagation of three optical fields in a 
triple mode optical fiber by applying the asymptotic reduction procedure [|2FI . Eqs.©-© has 
also been derived to describe the interaction between a quasi-resonance circularly polarized 
optical pulse and a long-wave electromagnetic spike [|29l . In the context of many-component 
magnon-phonon system, such a multicomponent YO system has also been proposed and its 
corresponding Hamiltonian formalism was studied [l30ll . Also, the authors in ref.[|28l have 
carried out Painleve analysis for Eqs.®-® and obtained the general bright N-soliton solution 
in the Gram determinant form. Soon after, they constructed an extensive set of exact periodic 
solutions in terms of Eame polynomials of order one and two lISTI . With the help of the KP- 
hierarchy reduction method, we have recently considered the dark soliton, mixed soliton and 
rational solutions for the coupled YO system in [[32l[33ll34]| . 

Integrable discretizations of integrable systems have received considerable attention 
recently (see, e.g., ifTTIl and references therein). So far, several approaches have been 
developed to construct integrable difference analogues of soliton equations. Ablowitz and 
Eadik proposed a method of integrable discretization by using Eax pairs [l35ll3^ . Based on the 
bilinear form, Hirota developed a method to construct integrable discrete analogues of soliton 
equations [|371[38l|39l. Date et al [|40ll4ni42ll43ll44ll have shown an effective way to discretize 
integrable equations via the transformation group theory, in which many of discretization 
versions of soliton equations can be derived by reduction. Suris [l45ll also developed a 
general Hamiltonian approach for integrable discretizations of soliton equations. In a series 
of recent research, the bilinear integrable discretization method has been widely applied to 
construct discrete schemes of the coupled Nonlinear Schrddinger equation [[T^fm . the (2-1-1)- 
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dimensional sinh-Gordon equation [|46ll , the two-dimensional Leznov lattiee equation [|47ll . 
the Camassa-Holm equation [l48l|43, the short pulse equation [l50ll and eoupled short pulse 
equation jSTl [52l|, the short-wave model of the Camassa-Holm equation |[53l . the redueed 
Ostrovsky equation [[54ll . the eoupled integrable dispersionless equation IfSSll . the integrable 
(2-i-l)-dimensional Zakharov equation [|5^ and so on. 

The goal of the present paper is to eonstruet an integrable semi-diserete analogue of 
the eoupled YO system by virtue of Hirota’s bilinear method, and derive both the bright and 
dark soliton (for the SW eomponents) solutions by the KP-hierarehy reduetion method. The 
remainder of the paper is organized as follows. In seetion 2, we present an integrable semi- 
diserete version of the eoupled YO system. In seetion 3 and 4, the bright and dark soliton 
solutions in terms of pfaffians of the semi-diserete eoupled YO system are eonstrueted based 
on two types of Baeklund transformation of the BKP hierarehy. Seetion 5 is eoneluded by 
some eomments and diseussions. 

2. Integrable semi-discrete coupled YO system 

Through the dependent variable transformation 

= L = 2(lnF)„, (6) 

F 

the ID eoupled YO equations ©-dS]) ean be east into the bilinear form 

(iA-D2)G(^)-F = 0, a^=1,2,---,M, (7) 

M 

D^DfF ■ F - 2cF ■ F = 2 G^^) , (8) 

A'=l 

where c is an integral eonstant and “ means eomplex eonjugate. The Hirota’s bilinear 
differential operator is defined by 

(|-T) 

By diseretizing the spaeial part of the above bilinear equations, 

^ + (9) 

F,^^iF„+i-F„), (10) 

one ean obtain 

iDrGi^^ ■ Fn - - IG^^Fn + =0, = 1,2, • • • ,M, (11) 

1 ^ f \ -t \ 

-AT^+i ■ F„ - cF„ ■ F„ = £ d'^GTcT . 

^ ^,=1 


( 12 ) 
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Furthermore, we require that the diseretized bilinear forms are invariant under the gauge 
transformation 

F„ F„ expiqon) , -)■ exp(<5ro«), 

then one gets the gauge invariant semi-diserete bilinear YO equations 

iAGi'^^-F,-^(G;,'JiF„_i-2Gi'^V„ + G;,'^iF„+i)=0, = 1,2,• • • ,M, (13) 


^AA+1 ■ a -cF„+i ■ a = £ ^(G,'f^iGi''^ + Gi''^G,'f^i). 

^I=l 


Let 


M _^n j _ 

— 2 


2 / 

^2 

-*• r? 


- 1 1 , L„ = 21nF„, 


(14) 


(15) 


the bilinear equations (fT3]) - (fT4l) are transformed into 




1 M 

hh+l.. -t.,) -2c = £ c<'''(sW.5W +sWsM,). 


(17) 


A'=l 


By using the relation ^ (L„+i — 2L„ -|-L„_i) = ln( 1 -|- yL„), we propose the following discrete 
system 

I c(^) ^ /■ z'c(^) I c(^) 

dt " £2 


-^(Ci+^ri-25r)-27^„(5r^^i+5:^^4)=0, = 1,2,• • • ,M, (18) 


pi=l 


dt 


(19) 


which converges to the coupled YO system dH)-© when £ —;■ 0. 

For simplicity, by taking £ = 1 and applying the gauge transformation S)^' —)■ 
exp(2iA)5i'“\ one can obtain the semi-discrete YO equations 






M 


or 


/y=l 

/ ^ -' yj \ 1 ^ ^ ^ \ 

lll(l + ~^) = 2 ^Ln+l — 2Ln+Ln-\), 


= (jM 1 + 1 ) (1 + iz,„), ^ = 1 , 2 ....,M. 


dt 

d 

dt 


M 


2^1n(l + y) = £ +S'i‘'ddl-S't‘-,)] 

In the subsequent two sections, we will consider general bright and dark soliton solutions 
for semi-discete coupled YO system (l23]) - (l24l) in details. For brevity, we refer to the soliton 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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solution including bright or dark soliton for the SW components and bright one for the LW 
component as the bright or dark soliton solution. 


3. Bright soliton solution for the semi-discete coupled YO system 


In this section, we will construct bright soliton solution for the semi-discete coupled YO 
system (l2^ - (l24l) . First, we briefly recall Backlund transformations of the semi-discrete BKP 
hierarchy by the following Lemma iTTOll . 

Lemma 3.1 The following bilinear equations 


Dtg^n^ ■ fn = gtilfn-i “ g^nljn+\ , 

■ fn = h^aljn-i - , 

-fn = g^nl\h\^^ -gn^h^nll, 

for /j = 1, • • •, M are satisfied by the pfafflans 

fn = Pf(«l, • • • ,<32Af,C2Af,'' ' ,Ci), 

g\f^ =pf(Jo,«i,---,«2Af,c2iv,---,ci,a^‘“)), 

= pf(Jo,6!l, • • ■,a2N:C2N, ' ' ' 
where the pfaffian elements are defined by 

pfiaj^at) = -^^—^(pjPkTexpi^j + ^k), 
pf{di,aj) = p"^^exp{^j), pf{aj,Ck) = 

pf(^i?/,c,■) =pf(J/,a(^)) =pf(J/,(3(^^) =pf(a/,a('")) = pf(a,-,p('")) =0, 


^ exp(Cf+nl' 


pf(c,,Q) = l e“+<> 

0 , 


7, VA/ J 7, 

^ 2N>j>N+l,N>k>V, 


otherwise. 


pf(cy,a(^)) = 
pf(c,-,pW) = 


exp(ri5^)), N>j>V, 


0 , 

0 , 


2A^> j >A^+1, 
N>j>U 


exp(^f), 2N>j>N+f 


with 




Here pj, ■’^jQ arbitrary constants. 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 
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Proof From the definition of the funetions /„, and we ean derive the following 
pfaffian’s rules: 

f„+i =pf(do,di,»), f„-i =pf(do,d^i,»), dtfn = pf(d^i,di,»), 

s!!ii =pf(di,»,a^^^), =pf(J_i,«,a^^)), dtg!^^ =pf(do,d^i,di,»,a^f'^), 

^ 1+1 =pf(d~h»,^^^^), =pf(do,d^i,di,»,^^^^), 

where 

pf(Jo,^i?i) =pf(^:?o,<i-i) = 1, pf(^:?-i,^:?i) = pf(a*^^\(3^-"^) =0, 

and (•) = (ai,---,a 2 N,C 2 N,---,ci)- 

Now the algebraie identities of Pfaffian 

pf( Jo, ^i?-i, Ji, )pf(«) = pf( Jo, J-i, •)pf( Ji, •, 

- pf( Jo, Ji, •)pf( J-i, ) + pf( Jo, •, a^‘“))pf( J_i, Ji, •), 

pf( Jo, J-i, Ji, •, )pf(«) = pf( Jo, J-i, •)pf( Ji, •, 

- pf( Jo, Ji, •)pf( J-i, •, p(^^) + pf( Jo, •, p^^^ )pf( J-i, Ji, •), 

and 

pf(Jo,Ji,«,a(^\p('"))pf(«) =pf(Jo,Ji,«)pf(«,a(^\p^'“)) 

-pf(Jo, •, a(^^)pf(Ji, •, p^-")) +pf(Ji, •, a^'“))pf(Jo, •, p^-"^), 
together with the above Pfaffian expressions of tau funetions give the bilinear equations (l25l) - 


Assume c = 0 in Eq. (l2T]) . whieh implies the bright-type soliton solution for semi-diseete 
eoupled YO system. Through the dependent variable transformation 

cip) _ ( j 

— A /• ’ \ r "> ^ 

Jn Jn 

Eqs. (I20l)-(l22l) are east into 

,(i“) ^ f — „(a') f ^ Ad) 


V fr 


n—\ 


DtgA •fn= 8n+lfn-l “ l/«+l, 

M 


Dtg\i^ •fn = -gi^i/„+l. 


Dtfn+l*fn = Y, - 8n^ g^n+l) ^ 

^=1 

for fj= 1,---,M. 

In order to carry out the dimension reduction, we define 

pf(a},4) = ^\ {PjPkT, 

P jPk A 


1), L = 21n/„, (31) 

(32) 

(33) 

(34) 
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pf(J/,a') = pf^,4) = Sj/t, 

pf(^i?/,c'■) =pf(J/,a(^)) =pf(J/,(3(^^) =pf(a'-,a('")) = pf(a'-,p('")) =0, 

g. 1 lN>j>N+l,N>k>V. 

Otherwise, 


pf(c;,4) = <; 

0 , 


pf(4,a(^)) = 


exp(ri(^4y, N>j>l- 


0 , 


2A^>7 >A^+1, 

pf(4,pW) = <| A^>7>1; 

1 exp(^5^^+^^-)> 2A^>7>A^+1- 

Then, the Pfaffians /„, and in Eqs.(l28l)-(l^ have alternative expressions in pfaffians 


fn — pf(4; ■ ■ ' 

= pf(do,a [, • • •, a'jjv, C2A,, • • •, 4, 

= pf(Jo, 4, • • •, 4iv> 4iv, • • •, 4, 

Therefore, under the reduction conditions. 


the following relation holds 


M 


Mfn.Sn\h^n^) = ^S^'^d^(^){fn,gn\h^n^), 

and thus one can get 

M 

‘n+l^ 


fl,/„+l ■/,. = E -slf’AS, 

/y=l 

Lastly, by taking the complex conjugate conditions 




= iV’ = PP ^S+i-;,o = for N>j>l, 


M 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


and y(^) are pure imaginary, then the function = gll‘\ and Eqs.(l25l). (l26l) and (1^ become 

Eqs.dH-dH. 

To summarize, we arrive at a Theorem below: 

Theorem 3.1 The bilinear equations (l3^-(l34l) are satisfied by the Pfaffians, 


fn = pf(ai, • • •, 

/ / 

4), 

(42) 

=pf(Jo,«i 

( ( 

r"j^2N^^2N 

,---,4,a(^4 

(43) 

=pf(Jo,«i 

( ( 
,Cl2N'=^2N 


(44) 
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where the pfaffian elements are defined by 

pf(a'-,4) = -^^—^{PjPkT, 
P jPk i 

pf(J,,a') = 


with 


pf(J,,c') =pf(J/,a(^)) =pf(J/,(3(^)) =pf(a'-,a('")) = =0, 

Pf(4’4) = 


2-u=i 


(rt I pM 




0 , 


2N>j>N+l,N>k> 1; 
otherwise, 


pf(c'-,a(^)) = 
pf(c',pW) = 


exp(ri;^'+y, N>j>V, 

0, 2A^>j>A^+1, 

0, N>j>U 

oxp{&^+^j), 2N>j>N+l. 


n 




Q 




forA^>y>l; 

2A^>7>A^+1, 

pj 


t.-(n _Lw 

^j-\Pj p.lP S' “S.O’ 


and Pj, and are eonstants satisfying p2N+i-j = Pj, ^2N+i-j,o = 

for > 7 > 1. 

In what follows, we will illustrate one and two bright soliton solutions for M = 2. 
One-soliton solution: By taking = 1 in (|4^ - (|4^ . we get the tau funetions for the 
one-soliton solution 


_ i^^UiAi+c(2)BiBi)pipipi-pi 

Jn- +2 


(PiPi-1)^ Pi+Pi 

=Aip"exp(^i), =Bipiexp(^i), 

where = (pi - = exp^j^^] and Bi = exp[p[^^]. 


— (piPi)"exp(^i+^i), (45) 

(46) 


In order to avoid the singularity, the eondition > 0 need to be 


2(P1+Pl) 

satisfied. Further, the above tau functions lead to the one-soliton solution as follows 
Sn'’ = y exp(i4/-0o)sech(^i^-h0o), 

= yexp(i4/-0o)sech(4i; + 0o), 


(47) 

(48) 


— 


(PiPi-i; 


sech2(^S + 0o), 


2pipi 


(49) 







10 


where + = n\n{ipi)+l,i and exp(2eo) = i . The 

quantities J^exp(—Go) and ■l^exp(—Go) represent the amplitudes of the bright solitons in the 
SW eomponents and S^n'^ respeetively. The real quantity denotes the amplitude 

of soliton in the LW eomponent. As an example, we illustrate one-soliton in Figure 1 for the 
nonlinearity eoeffieients = (1,-1). The parameters are ehosen asAi = l,5i = l+i 

and Pi = 1 + i. 




Figure 1. The profiles of evolutions of one-soliton solution (bright soliton for SW 
components). 


Two-soliton solution: By taking A = 2 in (I42l) - (l4^ . we get the tau functions for the 
two-soliton solution 

fn = l+CiiEiEi-\-Ci2E\E2 + C2lE2Ei-\-C22E2E2, 

+ \Pi2\\PnP22CnC2l-Pi2P2lCnC22)EiE2EiE2, (5G) 

=Ai£’i +A2£’2+^i2(Ai/’ijC2i — A2/2lCn)£’i£’2£i 

+ ^ 12 (^ 1 ^ 12^22 -MP22El2)ElE2E2, (51) 

gP =-Bi£’i +B2A2+^12(5l/’ijC2i — B2^2lQl)^1^2£l 

+ Pi2{B\P\2C22~ P'2.P22 Ei2)E\E2E2, (52) 


Pi -P2 p ^ Pj-Pk 
P1P2-1’ PjPk-l' 


with 
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Ej=p’]tx^{^j), Cji 


1 {c^'^'^AjAk + c^'^'^BjBk)pjPkPj-Pk 

2 {PjPk-^Y Pj+Pk 


where — {pj — y)t, Aj = exp[riyQ] and Bj = exp[riy^Q] for j = 1,2. Figure 2 shows two- 
soliton interaetion with the parameters as = ( 1 , —l),Ai = 2 ,A 2 = 5,Bi = 3,52 = 4, 

Pi = l+i andp 2 = -5 + fi- 




Figure 2. The profiles of evolutions of two-soliton solution (bright soliton for SW 
components). 


4. Dark soliton solution for the semi-discete coupled YO system 

In this section, we will consider dark soliton solution for the semi-discete coupled YO system 
(i23])- (l24l) . To this end, we need to introduce another set of Backlund transformations of the 
semi-discrete BKP hierarchy by the following Lemma IfTOll . 

Lemma 4.1 The following bilinear equations 

(A + OC^^^ — ' fn = ~ (53) 

~ <54) 

- a<^t + ^)/»+i ■ /, = iA)"' + Jjy S"'’aSi. <55) 
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for = 1, • • • ,M are satisfied by the Pfaffians 


f =t0-0 „(a') 

Jn '^n ’ <5^ "> ’ 


(56) 


where the tau function x 




is the Pfaffian x 


z(i).../(") 


= pf(ai,a 2 , • • • whose elements 


are defined by 


pf(a j,ak) = Cjk+ ^^ pf(Jo,a ;)pf( Jo, «/t), 

P 7 PZ" i 


M 


pf{di,aj)= p']+^Yl T 

^/=l \ ^ 








|/=1 
M 


■ p 


Pj 


7, 


exp(^;), 

1 — Pj 




• a(^) pj Pj — a(^) 

Proof From the definition of the tau functions, we have the following formulae of 
pfaffians: 


Wi 


= pf(Jo,Ji,»), =pf(Jo,J-i,»), 


= pf(Jo, Jl, d ^\, jj^^, •), 


/(1).../M + 1.../(M) _ 


= pf(Jo,j[^^«), 


^ ^n+l ~ 

1 /(i).../(/')+i.../W 

a(/') 


= pf(Ji,j|^\«), 

= pf(J_i,jj^\«), 


+ =pf(Jo,J-i,Ji,jj^^«), 

xr’-^"^-^-“'=pf(Jo,ji^J,.), 

1 


a(^) 


where 


pf(Jo,Ji) =pf(Jo,J-i) = 1, pf(J-i,Ji) =0, 
pf(ji^j\j|^^) =0, pf(Jo,j|^^) =pf(Jo,ji^i) = 1, 
pf(Ji,jj^^) = a(^), pf(Ji,ji^J) =pf(J_i,jS^^) = 


pf(j}^\a^) =p 


Pj 


a 




M 


\ 1 — a (^)p 


J v= 


n I 

V=1 V ^ 


Pj 


a 


(V) 


z(v) 


-aMp, 


exp(4), 


and (•) = (ai,---,a27v). 
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Now the algebraic identity of pfaffians together with the previous rules for tau functions 
in pfaffians 


pf( Jo, , d \, , •)pf(») = pf( Jo, J-i, •)pf( Ji, df ^, •) 

- pf( Jo, Ji, •)pf( J-i, jj'"^, •) + pf( Jo, j|^^, •)pf( J-i, Ji, •), 
gives the bilinear equations (l5^-(l54l) while the algebraic identities of pfaffians 

pf( Jo, Ji, ji^J, df ^, •)pf(«) = pf( Jo, Ji, •)pf( ji^J, df ^, •) 

- pf( Jo, ji^J, •)pf( Ji, d ^^, •) + pf( Jo, jj'"^, •)pf( Ji, d ^\, •), 
leads to the bilinear equation (l55]).D 

Next, we consider c 7 ^ 0 in (| 2 TI) . which implies the dark soliton solution of the semi¬ 
discrete YO system. By applying the dependent variable transformation, 

sr = (ia(^))”exp[(a(^) ST = (-ia(^))”exp[(a(^) 


Ln = 2 ' - 1^ , i = 21n/„ (57) 

with = 1, Eqs. (I20l)- (l22l) are cast into 

(A + ./„ = , (58) 

(A + ./„ = g^!!ljn+i , (59) 

A/«+l •fn-cfn+l»fn = ^ gl^hgn^ “ gi^^ gl^"} l) ^ ( 60 ) 

for fi= 1,---,M. 


Now we carry out the reductions to obtain bilinear equations (l58])-(l60l) from Eqs. ([53])- 
(1551) in Lemma 4.1. Eirst, by taking 

Cjk = ^ 2 N+l-j,k, j<k, (61) 

the tau functions can be rewritten as 

IN 

;=i 

where 


1 


pf(a'-,a^) — ^2N+i-j,k —^— x ./, 

J pf{do,aj)pf{do,a2N+i-j. 

Thus, if pj satisfies the constraint condition: 


PjPk-^ 


M 

fj=l 


i/j) ( _ 1 P2A'+1-; _ 1 -a(^V2Af-H- 


1 - a^dpj pj- aW 1 - a^dp2f^_^i_j p2N+i~j - 


1 

— Pj -|-P2A^+l-i “ 

Pj 


1 

P2A'+1-/ 


(63) 


( 64 ) 
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i.e., 


M / 1 


/y=l 


then we have 


aW 


{l-a^)pj){l-^^P2N+\-j) {\-ai^‘)p2N+\-j){l-^Pj) 


1 1 

—- ^ -. 

Pj PlN^l-j' 


^i=i 

Applying the above relation (1^ to Eq. (1551) . we have 


M 




M -i M 1 

\D, - £ = - E - ^4"'*^,). 


^=1 ^=1 
Furthermore, imposing the complex conjugate conditions 

-M 


(65) 


( 66 ) 


(67) 


V’ ^ 2 yv+i-; = > 7 > 1 , ( 68 ) 

and requiring being pure imaginary, I = one can have the function = gl(‘\ 
Consequently, Eqs.(l53l).(l54l) and (1671) become Eqs.(l58l)-(l60l). In summary, the following 
Theorem holds: 

Theorem 4.1 The bilinear equations (I58l) - (l60l) with |a(^) | = 1 are satisfied by 


(69) 


f -j.o- 0 Ap) _ ...0 AA 

Jn — ^ &ri — ’ 


where the tau function is the Pfaffian = pf(ai,a 2 , • • • ,« 2 A^), whose elements 

are defined by 

pf(<3;,«/t) = hN+i-j,k+ ^\ pf{do,aj)pf{do,ak), 




PjPk - 1 


M 


pf(Jy,a,)=pf^nll^ "" 

Pj 


A) 


l(p) 


aA) p . 


exp(^,), 


Here pj, ^jq are constants satisfying the following constraints. 

M AA 


■ Y. E{a^^\pj,p2N+i-j) = — H- - 

^tl 2 V J Pj P2N+l-j 

1 


with 
.X'YryM 


( 70 ) 


( ,Pj,P2N+l j) 
and p 2 iV+i -7 = Pj, 4 S+i„ 7 ,o = for A > 7 > 1. 














15 


In the following, we will illustrate one and two dark soliton solutions forM = 2. 
One-soliton solution: By taking = 1 in (l69l) . we get the tau funetions for the one 
soliton solution: 

. Pi -Pi 


/n — 1 + i 


-(;?ipi)"exp(^i+^i), 




(2) . , . Pi-Pl 

gn - 1+1- 


PlPl-l 

(pi, _,„ g, 

- tn-nr^Ti- 

pipi - 1 (1 -a^iVi)(l 

(pi-a(^))(pi g, 

MpiPiY 


PiPi - 1 (1 -a(2)pi)(i -a(2)pi) 

where = {pi — ■^)t + ^yo and pi is a eomplex eonstant satisfying 
( aM_a(A')^ f -- ^ ^ 


(71) 

(72) 

(73) 




A'=l 


1—a(‘“)pi)(l—a(^)pi) (1—a(i')pi)(l—a(i')pi^ 


1 1 

= —+—.(74) 
Pl Pl 


In order to avoid the singularity, the eondition > 0 need to be satisfied. Further, 

the above tau funetions lead to the one-soliton solution as follows 


c(l) _ ^ „i(n(pi + Slt+2fsin(pi) 

c(2) _ 1 i(n(p9 + fit+2fsin(p2) 

" - 2^ 

Ln = seeh^(^i^ + 9o), 

2pipi 


l+g24i 


) tanh(^i^ -h Go) 


e2^‘f’2)tanh(^i^-heo) 


where^'j =^'i^-fi^'i^ = nln(ipi)-h^i,exp(i(pi) =a(^^,exp(i(p 2 ) = a^^), exp(2eo) =i 


exp(2i(|)i) = 


(Pi-«(P)(pi-aP)) 


(l-a7)pi)(l-a7)pi 


A, exp(2i(^2) = ^ 


_ (pi-a+)(pi-a+) 

-a(2)pi)(l-a7)pi; 


(75) 

(76) 

(77) 

Pi-Pi 
pipi-i’ 

and the parameter pi is 

determined by Eq. (1741) . For the dark soliton in the SW eomponents and S^\ their 
intensities approaeh 1 as n —)■ ±oo, and the intensities of the eenter of the solitons read 
eos(|)i and eos(|) 2 . The real quantity denotes the amplitude of the soliton in 

the LW eomponent. We illustrate single dark soliton for the ehoiee of the nonlinearity 
eoeffieients (c*^^\c(^l) = (1,-2) in Figure 3. The parameters are ehosen as ^ + 

a(2) = 1 + ^21 = 1 - 0.7796i and ^i,o = 0. 

Two-soliton solution: By taking N = 2 in (l69l) . we get the tau funetions for the two- 
soliton solution: 


fn — i+Ei+E2 + 


{Pl -P2){Pl -Pl) 


gi ^ — 1 -\-^iE\ +A 2 E 2 + 

= 1 +EiE\ +B 2 E 2 + 


(Pl7’2-l)(PlP2-l) 

{Pl-P2){Pl -Pl) 


(P1P2-1)(P1P2-1) 
{Pl-P2){Pl-Pl) 


{piPl-i){piP2-l) 


EIE 2 , 

(78) 

Pl) 

2 

A\A2E\E2, 

(79) 

-1) 

Pl) 

1 

BIB 2 EIE 2 , , 

(80) 

-1) 
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Figure 3. The profiles of evolutions of one-soliton solution (dark soliton for SW components). 


where 


Ej = i (pjPj) exp(^,- + 1^), 

PjPj-l 

’’ ^ (1-a(2)j9y)(i-a(2)p^-) ’ 

and = {pj — y)t + and pj are eomplex eonstants satisfying 
_ y r_ I _^_ I _1 = 

2 V / (1—a(f')py) (1—a(f')py)(l— 


1 1 
— 

Pj Pj 


( 81 ) 


for j = 1,2. Figure 4 exhibits sueh a two-soliton interaetion with the parameters as 
(c(i),c(2)) = (i,-i), a(i) = ^ + ^i, a(2) = 1 - Pi = 1 -0.6714i, p2 = 0.5 - 1.5356i 
and ^'10 = ^^ 0 = 0. 


5. Conclusion and discussions 

In the present paper, we construet an integrable semi-diserete analogue of the coupled YO 
system by using Hirota’s bilinear method. Moreover, both the bright and dark soliton solutions 
in terms of pfaffians are derived based on the Backlund transformations of the semi-discrete 
BKP hierarchy. As far as we are aware of, it is the first time to propose an integrable semi¬ 
discrete YO system. There are several topics left for the study of the semi-discrete YO system, 
which are listed below: 
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Figure 4. The profiles of evolutions of two-soliton solution (dark soliton for SW 
components). 

(i) The Lax pair, the bi-Hamiltonian strueture and the eonservation laws of this newly 
proposed semi-diserete YO system are not elear and deserve to be explored. 

(ii) How to eonstruet the inverse seattering transform seheme for the semi-diserete YO 
system is also an interesting problem. 

(iii) It is of important to eonstruet the mixed-type soliton solutions with both the bright and 
dark soliton for the veetor semi-diserete YO system. 
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